Modern shell-model diagonalizations with realistic NN forces 
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The spectral and statistical properties of nuclei 46 V and 48 Cr are studied in the framework of 
nuclear shell model. A microscopical effective Hamiltonian derived from the CD-Bonn AW potential 
is employed. The calculations are preliminary results of a newly developed parallel shell-model 
code. Calculations with the realistic NN interaction reproduce well the ground-state rotational 
band and backbending phenomenon in 48 Cr. A sudden alignment of all particles is seen above the 
backbend. The nearest-neighbour spacing distribution in 46 V is presented and compared with the 
predictions of random-matrix theory. Significant deviation from the single Gaussian-orthogonal- 
ensemble distribution is obtained, indicating that the isospin-symmetry breaking is not so strong in 
the nucleus. 

PACS numbers: 21.30.Fe, 21.60.Cs, 24.60.Lz, 21.10.-k 



The nuclear shell model has been considered as the 
most fundamental tool for the investigation of nuclear 
structures and electromagnetic and weak transitions. It 
provides the microscopic basis to study the nature of 
nuclear states, magnetic and quadrupole moments, (3 
and 7 transitions and nuclear reactions. The simpli- 
est version of the shell model, the so-called indepen- 
dent particle model, has been proposed more than fifty 
years ago and gives successful descriptions of the nuclear 
magic numbers and ground-state angular momenta with 
the introduction of strong spin-orbit coupling interaction 
[l[ • Soon, the community realized the importance of in- 
troducing configuration-mixing in explaining the spec- 
troscopic properties of nuclei (i.e., 19 F) with singular 
ground-state spins 0, [|[ . The configuration- mixing con- 
cept has been considered as the essence of modern shell 
model calculations 0, H . 

Now the shell model has been well understood in the 
language of perturbation theory. The basis to study 
the nuclear many-fermion systems is a non-relativistic 
Schrodinger equation which has infinite freedom and in- 
solvable. As done in atomic many-body studies, the 
starting point is to reduce the Schrodinger equation to 
finite freedom with the help of Feshbach operators @, Q • 
The Schrodinger equation can be written in the form 
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where H is the original Hamiltonian, E the eigen ener- 
gies and P and Q projection operators defining the model 
space and the excluded orthogonal space. From Eq. (JTJ) 
we can derive an model-space dependent effective Hamil- 
tonian which is given by 0] 



H cS = P{H + H 
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The original eigen energies and wave functions can be 
constructed by diagonalizing this equivalent Hamilto- 
nian. The effective Hamiltonian in Eq. ([2]) is energy 
dependent and can be solved self-consistently Q or per- 
turbatively Q. In the latter case the effective Hamilto- 
nian can be written as 



H eS = PH P + PVflP, 



(3) 
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where Hq is the unperturbed Hamiltonian, Q the wave 
function operator and V the nuclear force. The last term 
is referred to as the effective residual interaction with 
v c g — PVVlP. In the shell model context, the single par- 
ticle Hamiltonian is evaluated from the experimental ob- 
servations of the single particle outside the assumed core. 
The effective interaction v c g is usually treated perturba- 
tively up to two-body level and expressed as two-body 
matrix elements in harmonic oscillator (HO) basis. The 
shell structure of nuclear systems recognized by Mayer et 
al. [l[ makes it possible to constraint calculations in one 
or several major shells defined by the projection operator 
P which can embed the dominant component of the the 
original wave function. 

Although calculations can be restricted in a single ma- 
jor shell, the corresponding Hilbert space dimension is 
still very large. The fp shell (with M-scheme dimension 
around 10 9 ) presents the state-of-the-art shell-model per- 
formance in the past decade. Diagonalizations for nuclei 
in the intermediate fp model space (e.g., 56 Ni), however, 
are still painstaking [lfj . Hence, the development and op- 
timization of shell-model diagonalization algorithm has 
been a longstanding work d M, El El El El El OS 
El Ell • The first implementation of efficient shell-model 
diagonalization is the Oak-Ridge Rochester code [ll| us- 
ing the j-j coupled scheme Another revolutionary 
development attributes to the Glasgow group |12]. The 
algorithm implemented by Whitehead et al. [12| intro- 
duced the M-scheme representation and the Lanczos di- 
agonalization method. The M-scheme is compatible with 
basic bit operations in modern computer and avoids the 
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FIG. 1: Schematic picture for shell-model wave function as 
an expansion of partitions. The lines indicate single-particle 
orbits in the HO basis. 



time-consuming calculations of coefficients of fractional 
parentage (CFP) in j-j scheme. The Lanczos algorithm 
uses an iterative procedure and is suitable for giant ma- 
trix diagonalizations [!, [3. The powerful shell -model 
code ANTOINE of the Strasbourg-Madrid group suc- 
ceeded the idea of the Glasgow algorithm [3| ■ 

In M-scheme, only J z and T z are good quantum num- 
bers, leading to a maximal dimension of the generated 
Hamiltonian matrix. The famous OXBASH code and 
its MSU version [l3| follows a hybrid algorithm between 
the M-scheme and the j-j coupled scheme. The sym- 
metry of J (T) is restored through a projection proce- 
dure. The dimension of the matrix can be significantly 
reduced through this symmetry restoration. The algo- 
rithm implemented in OXBASH makes it very efficient 
to calculate many converged vectors for a given angu- 
lar momentum. However, the projection process is very 
time-consuming and has precision problems. 

In this paper we report on preliminary results calcu- 
lated with a newly developed parallel shell-model code. It 
follows the hybrid algorithm of OXBASH and embeds a 
projection procedure to generate basis with good J quan- 
tum number. The program is written in Fortran 90 for- 
mat which is more efficient and is more suitable for high 
performance computing than the outmoded FORTRAN 
77 standard [l9(. In Fortran 90 language, the precision of 
data objects can be easily denned, while super-large ma- 
trices and arrays can be handled. The code is implanted 
on the 64-bit Beowulf cluster of PKU computer center. 
Parallel algorithms are realized with the help of Message 
Passing Interface (MPI) HJ. 

The starting point of the calculation is to generate 
a set of bases with good magnetic quantum number in 
M-scheme. Neutrons and protons are blocked into two 
spaces. The bases are classified according to the distri- 
butions of particles in the single-particle orbits (referred 
to as partitions). Schematic figure for shell model wave- 
function as an expansion of partitions is plotted in Fig. 
[TJ Since the angular-momentum projection operator can 
only change the magnetic quantum number, projected 
vectors in different partitions are orthogonal, ensuring 
that projections can be done for each partition separately. 
In Fig. [5] we plotted M-scheme dimensions, j-j scheme 
dimensions and number of partitions for + states in fp- 
shell N = Z nuclei as a function of neutron (proton) 
numbers. 

Vectors with good angular momentum are expanded 
in M-scheme bases for each partition, which can be given 
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FIG. 2: M-scheme and j — j coupled scheme dimensions for 
+ states in fp shell N=Z nuclei. Stars denote the numbers 
of partitions for each state. 



as 

M) = J2 M "nP J \a m ), (4) 

where P J is the projection operator and \a) a set of 
specially chosen M-scheme bases. M. is a lower trian- 
gle matrix. Vectors with good angular momentum are 
the orthonormal and satisfy, 

<¥f|tf/> = (*/| ^2 M jm P J \a m ) 

= 5ij. (5) 

Elements of M. can be obtained by inverting the matrix 
(^f\a m ). Since the projection operator satisfies, P J ■ 
P J = P J , the n-th vector can be obtained through, 

\Q n ) = M-^o.K) Y^Ma^), (6) 

i < n j < i 

and 

\* J n ) = P J \Q n )((Q n \P J \Qn))- 1/2 - (7) 

In the meanwhile, matrix elements in the n-th row of 
the M. is calculated. This projection algorithm is very 
efficient since calculating operations and the number of 
coefficients needed in the calculation are comparatively 
small and can be stored in the computer memory. Pro- 
jections for various partitions can be done independently 
and can run parallelly with very high scalability. 

The projection procedure shown above was suggested 
but not used in OXBASH 13] since it does not con- 
serve t z values and is unsuitable for projections in isospin 
space. The advantages of isospin space lie in that only 
the smaller set of (J,T) vectors satisfying J = Jm and 
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FIG. 3: Relative energies of first T = 1 states (J n = 0+) and 
first T — states for odd-odd N = Z nuclei (T z = 0) in the 
sd and O/7/2 shells. Experimental data are taken from Ref. 
231. Calculations for sd and O/7/2 shells are done with the W 
21 and CD-Bonn [U, KB3 Q and FPD6 interactions, 
respectively. 
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FIG. 4: Calculated and experimental yrast band in 48 Cr. 



conserved in the Hamiltonian, the two-body matrix ele- 
ments can be further simplified by adding isospin quan- 
tum numbers to the proton-proton and neutron-neutron 
(T = 1) and proton- neutron (T = and 1) interactions, 
which are restricted by J + T=odd for a = (3 and/or 
7 = S. The Hamiltonian matrix elements are given as, 

^i\H\^)=Y,M 3n ^{\H\a). (9) 



TABLE I: Examples of states in even-even nuclei that have 
T > \T Z \. Calculations are done using the W interaction [2lj . 



Nuclei 


r 


Cal. (MeV) 


Exp. (MeV) 


20 Ne 


It 


11.20 


11.26 


24 Mg 


1:1 


9.99 


9.97 


28 Si 


1 3 + 


10.81 


10.72 


32g 




7.06 


7.00 


36 Ar 


2.1 


6.65 


6.61 


36 Ar 


3t 


7.46 


7.34 


36 Ar 


1 3 + 


8.19 


8.13 



T = T z are projected. The number of non-zero Hamil- 
tonian matrix elements is also smaller. The restriction 
of J = Jm is valid since the Hamiltonian is invariant 
under rotations in angular momentum space. However, 
nuclear states with T > T z do exist and can exist at very 
low excitation energies. Examples in odd-odd and even- 
even nuclei are shown in Fig. [3] and Table HI respectively. 
To give the full level structure, besides the generation of 
the smaller (J, T) set with T = T z , projections in isospin 
space have to be imposed to the maximum of the T value. 

Since the nuclear force V is a scalar product and con- 
serves spin and parity, the effective interaction can be 
simply expressed as 

Wcff = E E (afcJ\V\j6;J), (8) 

J a<P;~/<8 

where a denotes the single-particle orbit. If isospin is 



In practical calculations, the matrix (^f\H\a) are calcu- 
lated first and store in memory. The total matrix ele- 
ments are generated by scanning and timing the M. ma- 
trix. The calculation of Hamiltonian matrix elements can 
be separated by partition and done parallelly. However, 
the efficiency is not optimal when dimensions are not very 
giant. Another very efficient way is to separate the above 
two steps into two or several computing processes. 

The Hamiltonian matrix is diagonalized using the 
Lanczos procedure The most time-consuming 

part is the iterative generating of Lanczos vectors, 

p nH )=H\0 n ), (10) 

where |/?o) is a specially chosen basis bxj-j scheme. In the 
present code, the Hamiltonian matrix elements are scat- 
tered into available processes using the MPI paradigm 
[20I ] . in which the matrix-array production can be done 
parallelly and be reduced to give out the new Lanczos 
vectors. 

The maximum capacity of the computer code has not 
been tested, which can be restricted by the Input/Output 
operation rate of the computing system. In this work we 
present calculations for nuclei 46 V and 48 Cr. The largest 
j-j scheme dimension is about 2.5 x 10 5 , for which the 
scalability of the diagonalization using up to 30 proces- 
sors is still optimal. 

In the past decade, the medium-mass nucleus 48 Cr, 
which has a half-filled O/7/2 shell of protons and neu- 
trons, has been under intensive study (see e.g. Ref. [27j 
and reference therein). A well-deformed ground state 
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FIG. 5: The percentage of pure O/7/2 configuration in the 
wave functions of yrast states in 48 Cr. 



appears in 48 Cr, together with many macroscopic phe- 
nomena, like rotation, backbend and triaxial deforma- 
tion [27], [28| . Previous shell-model study by Caurier et 
al. using a modified Kuo-Brown interaction reproduced 
the yrast band and the backbending phenomenon [28| . 

In the present work, We used an isospin-nonconserving 
effective Hamiltonian to study the structure of 48 Cr in 
full fp shell. The effective interaction is derived micro- 
scopically from a high-precision charge-dependent Bonn 
potential [261 ] . Single-particle energies are taken from ex- 
perimental observations of odd neutron and odd proton 
with respect to the core of 40 Ca, in which the Thomas- 
Ehrman shift of proton Ip orbits can also be taken into 
account approximately. Excellent agreements with ex- 
periments is seen till the band terminated state with 
/ = 16 + . Calculations characterize the backbend ap- 
peared at / = 12 + . Detailed analyse of wave functions 
shows that, for all states in Fig. [4] the average numbers 
of particles in 0/7/2 orbit (or O/7/2 orbital occupancy) 
are nearly constant (around six). Above the backbend, 
particles are suddenly aligned, which manifest itself in 
that pure 0/7/2 configuration significatcly increases for 
states with / > 12. The dispersion of configurations at 
low spins contributes to the deformation. 

The shell model provides the most promising founda- 
tion to investigate the complexity and chaotic behavior 
in nuclear many-body systems, which is the best test- 
ground for the study of quantum chaos. Comprehensive 
reviews on the topic can be find, e.g., in Refs. [29L l3fj|. 
Now it is believed that the level fluctuation properties 
(i.e., the nearest-neighbour spacing distribution (NNS)) 
of quantum systems with time-reversal symmetry whose 
classical analogs are chaotic can well described by the 
Gaussian orthogonal ensemble (GOE) of random-matrix 
theory (RMT), whereas quantum analogs of classically 
integrable systems follow the Poisson statistics. Partic- 
ularly interesting is the statistical properties of systems 
with one additional approximate symmetry. If the sym- 
metry is exact, the NNS distribution follows the super- 
position of two independent GOE (2-GOE) [29]. The 
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FIG. 6: The NNS distribution for low-lying levels in 46 V with 

J < : 




FIG. 7: Same as Fig. [6] but only for states with even spins. 



complete absence of the symmetry in the distribution, 
however, is expected even with a small breaking of the 
symmetry [3l|. A vivid example has been shown in the 
spectral statistics of acoustic resonances in rectangular 
quartz blocks, in which a transition from 2-GOE to GOE 
is seen with the gradual breaking of the point-group sym- 
metry in the system [33I ] . 

Systematical works have been done to study the 
chaotic behaviors of sd-shell nuclei, both experimentally 
and theoretically [30l| . The chaotic behavior in sd shell 
nuclei has been well described by shell model calculations 
using the empirical 



W interaction [2l|. However, con- 
troversy still exists concerning the role played by isospin 
symmetry breaking on the statical properties of NNS [3l[ . 
In nucleus 26 Al, a system with 5p — 5n above the core, the 
NNS distribution for low-lying states generated by em- 
pirical isospin-nonconserving Hamiltonian shows a strong 
isospin dependence, indicating that the isospin symme- 
try breaking is small, whereas a strong isospin symmetry 
breaking effect is expected from the analyse of available 
observations [3l[. In the sd shell, strong Coulomb in- 
fluence exists due to the large Thomas-Ehrman shift of 
proton lsi/2 orbit. This can lead to large uncertainties 
[3lj . One of the notable feature of the fp shell is that 
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the isospin symmetry breaking is relatively smaller [2(| . 
It would be very helpful to see the statistical behavior 
in the shell and the isospin symmetry breaking influence. 
Deviation from the GOE distribution has been noticed 
in the NNS distribution for the low-lying states of Ca 
isotopes [32|] ■ 

We evaluate the isospin symmetry breaking effect 
in the NNS distribution of 46 V using the isospin- 
nonconserving force [2(| discussed above. The lowest 150 
positive-parity states for each spin between 1 = and 16 
(O/7/2 band termination is at / = 15) were calculated. 
The 2550 states can be generated on a single processor 
within an hour. The energy levels for each spin are un- 
folded through [32j 

N(E) = [ p(E')dE' + N , (11) 
Jo 

where p(E) is the mean level density which is assumed 
to be a exponential form, 

p(E) = l/T-exp[(E-E )/T}. (12) 

The distribution of NNS, P(s), is obtained by accumu- 
lating the number of spacings of Sj = N(Ei + i) — N(E{) 
within (s, s + As). Results are plotted in Fig. [S] and [71 

Our calculations are carried out in proton/neutron 
space and is irrespective of isospin. The isospin sym- 
metry breaking in the effective Hamiltonian has three 
origins: the Coulomb shifts of proton lp single-particle 
orbits (~ 200 keV) , the charge dependence in strong force 
and two-proton Coulomb matrix elements. If isospin 
symmetry is exact, the distribution should follow the 2- 
GOE distribution, while a transition to the single GOE 



distribution is expected with the gradual breaking of 
isospin symmetry. Our calculations show that, however, 
the deviation from the 2-GOE distribution is relatively 
negligible. Other statistical properties, like the Brody 
parameter and the spectral rigidity, can provide more 
deep insight on this problem. 

In conclusion, a new high-performance parallel shell- 
model code has been developed. With the code we cal- 
culated the spectral and statistical properties of nuclei 
46 V and 48 Cr with a microscopical effective Hamiltonian 
derived from the CD-Bonn NN potential. Calculations 
with the realistic NN interaction reproduce well the yrast 
band and backbending phenomenon in 48 Cr. A sudden 
increasement of pure O/7/2 configuration is seen above 
the backbend. In the meanwhile, the average numbers 
of particles in the yrast states remains nearly constant. 
The NNS distribution for low-lying states in 46 V is cal- 
culated and compared with the predictions of random- 
matrix theory. The isospin-nonconserving Hamiltonian 
is used to evaluate the isospin symmetry breaking effect 
on the NNS distribution. Significant deviation from the 
GOE distribution is obtained, indicating that the isospin- 
symmetry breaking is not so strong in the nucleus. 
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